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ABSTRACT 

We calculate the neutrino pair annihilation rate into electron pairs between 
two neutron stars in a binary system. We present a closed formula for the energy 
deposition rate at any point between the stars, where each neutrino of a pair 
derives from each star, and compare this result to that where all neutrinos derive 
from a single neutron star. An approximate generalization of this formula is given 
to include the relativistic effects of gravity. We find that this inter-star neutrino 
annihilation is a significant contributor to the energy deposition between heated 
neutron star binaries. In particular, for two neutron stars near their last stable 
orbit, inter-star neutrino annihilation energy deposition is almost equal to that 
of single star energy deposition. 

Subject headings: gamma rays: bursts — gamma rays: theory 

1. Introduction 

Neutrino emission from, and interactions above neutron stars has long been understood 
to be an important process in supernovae (e.g. Wilson & Mayle 1993). One interaction of 
particular interest is neutrino annihilation: v+V — > e + +e~ because it converts neutrinos into 
electrons, which, due to their vastly larger cross- sect ion, are an important source of energy 
for supernovae and possibly gamma-ray bursts. Because of the strong gravity of neutron 
stars, relativistic effects must be included when calculating neutrino processes. Rates of the 
vv — > e + e~ interaction above a neutron star were originally calculated without gravity by 
Cooperstein et al. (1987) and Goodman et al. (1987). In Salmonson & Wilson (1999) these 
neutrino annihilation rates were calculated with gravity and it was found that general rela- 
tivistic effects significantly augmented the efficiency; by up to a factor of four for supernovae 
and up to a factor of thirty for lone neutron stars at their critical density. 
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Another situation where neutrino emission from neutron stars is important is in neu- 
tron star binary (NSB) systems. Binary neutron star mergers have long been a possible 
progenitor for gamma-ray bursts (GRBs) (e.g. Eichler et al. 1989; Meszaros & Rees 1992). 
In particular, current trends would suggest that NSBs may be the progenitor of the "short" 
GRBs, characterized by timescales < 1 second, which are thought to be distinct from the 
class of "long" GRBs (e.g. Fryer et al. 1999). As such, we are moved to study mechanisms 
for energy extraction. 

Due to the extreme high densities inside the neutron stars, neutrinos, with their small 
cross-section, are a necessary component of the energetics of supernovae and compact binary 
merger scenarios. In particular, only a small proportion (~ 1%) of the total energy (which 
could be as high as the binding energy of a solar mass neutron star) in neutrinos need be 
transferred to conventional matter in order to produce large, observable, energetic events. 
Such a mechanism has been suggested as an energy source for GRBs by Salmonson et al. 
(2001); Wilson et al. (1997); Meszaros & Rees (1992); Eichler et al. (1989). 

In this paper we study the inter-star neutrino annihilation; where each neutrino in 
an annihilation pair derives from each star in the binary. This is distinct from the single 
star neutrino annihilation discussed in Salmonson & Wilson (1999). Both processes will 
occur in a binary star system and will sum to give the total energy depostion from neutrino 
annihilation. We will compare and contrast the contributions of these two processes. We find 
that the inter-star component is a very significant contributor to the total energy deposition. 
In Section 2 we derive the results without general relativistic effects. In Section 3 we make 
some estimates of the augmentation of annihilation due to general relativity and compare 
with the relativistic augmentation derived in Salmonson & Wilson (1999) for the single star 
case. 



Starting from references in Salmonson & Wilson (1999), the energy deposition via 
neutrino-antineutrino annihilation into electron-positron pairs per unit time per unit vol- 
ume is 



where f v and f v are the respective neutrino and anti-neutrino number densities in phase 
space, e v and e v are their energies, 6(f) is the integral over incident neutrino angles, G\ = 



2. 



Interstar Neutrino Annihilation 




(1) 
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5.29 x 10" 44 cm 2 MeV" 2 and 



D = 1 ± 4 sin 2 O w + 8 sin 4 9 



w 



(2) 



where the + sign is for v e V e pairs, while — is for v^y^ and v T V T pairs. We assume black- 
body-like neutrino emission; a good assumption for thermal neutrinos diffusing out of a hot 
neutron star (Wilson et al. 1996). Neglecting redshift effects on the temperature for the 
moment, eqn. (1) reduces to 

q = 7.97 x 10 21 D 0(f) ergs cm" 3 sec" 1 (3) 

\ lMev / 

where D = 1.23 for v e V e pairs and otherwise D = 0.814. The integration over incident 
neutrino angles is 

0(f) = ff(l-n u - SljjfdVLvdVLjj . (4) 



We define the neutrino trajectories, in (x,y,z), from star 1 (on the left) as 

fii = (cos#i, sin #i cos 0^, sin 9[ sin (5) 

and for star 2 (on the right) as (note that x — > —x do to reflection symmetry) 

f2 2 = ( — cos# 2 , sin# 2 cos0 2 , sin0 2 sin0 2 ) (6) 

where (9' 1 ,4>' 1 ) and {9' 2 ,<f) 2 ) are the direction angles for the neutrino trajectories in spherical 
coordinates from stars 1 and 2 respectively, so that {9[ = 0, $ = 0) points to the center of 
star i — 1,2. Since we wish to evaluate the integral of eqn. (4) at any point in space around 
the stars, we rotate the vector by an angle Q 1 

$7i = (cos 9[ cos ©i — sin 9[ cos (/)[ sin 6i, sin 9[ cos cos Q 1 + cos 9[ sin 1; sin 9[ sin (/>[) . 

(7) 

A rotation of fl 2 by 6 2 gives (note that f2 2 will have the same form as fli except x — > —x) 

f2 2 = (— cos 9' 2 cos 9 2 + sin 9' 2 cos (f> 2 sin 2 , sin 9' 2 cos 2 cos 2 + cos 9' 2 sin 6 2 , sin 9' 2 sin 2 ) . 

(8) 



Thus, noting (1 - Cl l ■ Sl 2 f = 1 - 2fti • ft 2 + (^i ■ ^) 2 and dO* = sm9' i d9' i d(j)' i where 
i = 1,2, eqn. (4) can be integrated and becomes 
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©inter(0i, 2 , 0i, 2 ) = (vr 2 (516 + 8 cos[0i] + 4 cos[20i] + 4 cos[0i - 20 2 ] + 8 cos[0i - 2 ] + 2 cos[2(0i - 2 )] 

+ 4 cos[20i - 2 ] + 8 cos[0 2 ] + 4 cos[20 2 ] + 8 cos[0i + 9 2 ] + 2 cos[2(0i + 2 )] + 4 cos[20i + 2 ] 

+ 4 cos[^ + 20 2 ] + 12 cos[0i - 20! - 6 2 - 20 2 ] + 6 cos[20i - 2@i - 9 2 - 20 2 ] 

+ 12 cos[0i - 2©x + 2 - 20 2 ] + 6 cos[2#i - 2@x + 9 2 - 20 2 ] + 96 cos[0i - ©i - 2 ] 

+ 6 cos[2(0i - 0i - 2 )] + 48 cos[0i - ©i - 2 - 2 ] + 3 cos[2(0i - ©i - 9 2 - © 2 )] 

+ 48 cos[0i - ©i + 9 2 - @ 2 ] + 3 cos[2(#! -Q l + 9 2 - © 2 )] + 192 cos[©i + © 2 ] + 12 cos[2(©! + © 2 )] 

+ 96 cos[0i + ©i + 2 ] + 6 cos[2(0i + ©i + © 2 )] + 96 cosf©! - 9 2 + 2 ] + 6 cos[2(6i - 9 2 + © 2 )] 

+ 48 cos[0i + ©i - 9 2 + 2 ] + 3 cos[2(0i + ©i - 9 2 + © 2 )] + 96 cosf©! + 9 2 + © 2 ] 

+ 6 008(2(0! + 9 2 + © 2 )] + 48 cos[#i + ©i + 2 + 2 ] + 3 cos[2(0i + ©x + 9 2 + © 2 )] 

+ 12 cos[2©i - 9 2 + 2@ 2 ] + 12 cosf^ + 2@x - 9 2 + 20 2 ] + 6 cos[20x + 20 : - 9 2 + 20 2 ] 

+ 12 cos[2©i + 9 2 + 20 2 ] + 12 cos[0i + 26i + 2 + 20 2 ] + 6 cos[20i + 2©i + 9 2 + 20 2 ] 

+ 12 cos[0i - 2(0x + 2 )] + 12 cos[0i + 2(©i + 2 )] + 6 cos[#i - 2(0 X - 9 2 + 2 )] 

+ 6 cos[0! + 2(0! - 9 2 + 2 )] + 6 cos[#i - 2(0 X + 9 2 + 2 )] 

+ 6cos[0i + 2(0! + 2 + 2 )]) sin 2 [0i/2] sin 2 [0 2 /2])/24 . 



An illustration of the geometry is shown in Figure 1. Specifically, 0« is the apparent ra- 
dial angular size of star i and ©« is the angle star i is observed at with respect to the 
axis connecting the star centers. It is useful to convert eqn. (9) to cartesian coordinates 
(0(0i, 2 , ©i, 2 ) — > <d(x, y)) via the following substitutions. If two stars of mass M and 
radius R are separated by distance d, so the half-separation is d 2 = d/2 and 



n,2 = V (d 2 ± x) 2 + y 2 , 



(10) 



then 




(11) 




(12) 
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where the + and — signs correspond to stars to the left (star 1) and right (star 2) of the axis 
of symmetry respectively. Figure 2 shows a contour plot of this integration over inter-star 
incident neutrino angles given by eqns. (9, 11, 12). As expected, vv annihilation takes place 
primarily between the stars, where the cross-section and flux are highest. 

For comparison, the result for a single neutron star is (Salmonson & Wilson 1999) 

2tt 2 

0(r)sin g ic = — (1 - x)\x 2 + 4a; + 5) (13) 

where r is the distance from the neutron star center and, for a star of mass M and radius 
R, we define 



-^-(tJttw- (14) 

This equation (14) takes into account the effects of general relativity in Schwarzchild coor- 
dinates. In particular, the gravitational bending of neutrino trajectories (see Section 3). For 
the Newtonian result, take M = 0. In Figure 3 is shown a comparison between 0i n t er (eqn. 
9) and S ingic (eqn. 13) along the axis between the stars. We see that annihilation of vv from 
a single star is dominant near the stellar surfaces, while inter-star annihilation is important 
between the stars. 

One can compare the rate of inter-star deposition to single star deposition. A key 
conclusion of this paper is that the inter-star annihilation is a significant component of the 
total energy deposition. In Figure 4 we see a ratio of total energy deposited via inter-star 
annihilation versus that for a single star over a range of stellar separations. The total energy 
deposited is an integral of eqn. (3) over the interstellar volume V 



Q 



Newt 



J qdV oc J dV . (15) 



We see that in this Newtonian limit Qmter/Qsmgie ~ 80% for close separations (i.e. c?2 < 2/2) 
where Qinter oc / int cr dV using eqn. (9) and Qsingic oc / O sing i e dV using eqn. (13). 



3. General Relativity 

Analytical calculation of the general relativistic effects on inter-star neutrino annihi- 
lation is extremely complex. Instead, we make some crude approximations to gain some 
insight into the effect of gravity. We neglect possible rotation effects. It is most convenient 
and intuitive to use isotropic coordinates 

ds 2 = -a 2 dt 2 + <j)\dr 2 + rW) (16) 
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where, for one star (Misner et al. 1973) 



1 _ M 
a = - 2r 

i , M 

I (17) 

and the isotropic distance, r, to the center of the star is related to the Schwarzchild distance 
r by 

r = r^l + ^0 =F0(r) 2 . (18) 

In the two star environment we define 

. . MM 

a(r+,r_) = 1 - — 

M M 1 yj 

where r + (r_) indicates the distance to the center of the star on the left (right). This choice 
of metric is found to be an adequate approximation, to a few percent, of the real metric as 
numerically calculated by Wilson et al. (1996). The value of these metric components at 
the stellar surface is not at constant potential (because we neglect tidal distortions which 
equalize the surface potential), thus we define the potentials at the surface to be 

a s = a(R, d s ) , <p s = <j>(R, d s ) (20) 

with 




l)±x) +y 2 (21) 

using eqn. (10) and where one chooses the upper sign if x > and vice versa. This corre- 
sponds to the potential at the surface of the star at the point for which the line to point 
(x, y) is perpendicular to the surface. From here on we simplify the notation by suppressing 
the metric arguments (e.g. a = a(r + ,r_)). 

For the single star case, in these coordinates, eqn. (14) becomes 

(22) 
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where R is the isotropic stellar radius. For compact neutron stars near collapse, we will take 
R = 7.5 km (Wilson et al. 1996). 

For the inter-star case, the apparent radial angular size, 9, of a star will be augmented 
due to the gravitational bending of neutrino geodesies by the mass of the star. We include 
this effect by modifying eqn. (11) by identifying R — ■> RcjfJ 2 /at s and r — > r<f) 3 ^ 2 /a from 
eqn. (22). These substitutions take into account the "lensing" effect due to the gravitational 
bending of the neutrino geodesies (Salmonson & Wilson 1999) and give 

■ {Uj ffe)) ' (23) 

In addition, there will be a redshift (or blueshift) of the neutrino energies. Specifically, 
assuming the neutrinos emerge from the neutrinosphere of each star with surface temperature 
T s , then the temperature at a position (x,y), using eqn. (10), is 

T(r + ,r_) = ^T s . (24) 
a 

Then the energy deposited at this point is (eqn. 3) 

q = 7.97 x 10 21 D ( -^rrrl 0(f) ergs cnT 3 sec' 1 oc ( — t}\ 6(f) . (25) 
\ lMev / \ a J 



Following Salmonson & Wilson (1999) we define the total energy deposition 

Qgk = j ' qdV = 7.97 x 10 21 ,D [^^j J ^^ dx " er S s sec_1 ( 26 ) 

where the metric is given by eqn. (19) and the inter-star angle factor is calculated using 
eqns. (9, 12, 23) and the single star angle factor is given by eqns. (13, 22). A comparison of 
these two rates for several stellar masses and a range of separations is shown in Figure 5. One 
can see that as the stars approach their last stable orbit, d 2 ~ 10 — 12 km (Mathews & Wilson 
2000), the inter-star contribution to neutrino energy deposition, Qinter becomes comparable 
to the single star contribution Qsingic- Thus the inter-star energy deposition is significant for 
hot neutron stars in close binaries. In Figure 6 is shown the augmentation of both inter- 
star, Qinter, and single star, Qsmgie, energy depositions for a typical half-separation, d<i = 15 
km, and a range of masses. These are compared with the energy deposition for an isolated 
neutron star, Q one star , (Salmonson & Wilson 1999) here presented in isotropic coordinates. 
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4. Discussion 

Herein we have calculated the inter-star neutrino annihilation rate for close neutron stars 
of equal mass in a binary system. A generalization of interest would be the consideration 
of systems with stars of unequal mass. While a detailed calculation of such is difficult and 
beyond the capability of the analytical results presented here, we can argue that, at least 
for small departures from mass equality, the inter-star annihilation rate will not change 
much. In the stellar compression effect discussed by (Mathews & Wilson 2000), the larger of 
the two unequal mass stars may be thought to receive larger compression due to its excess 
mass. From the results given in tables III and IV of Mathews & Wilson (2000) is found 
that the compression for a given four-velocity is quadratic in mass. However, the smaller 
mass star will move faster, varying inversely with mass, and compression is shown to scale 
quadratically with four- velocity (Mathews & Wilson 2000). Therefore we expect that there 
will be a balancing effect keeping the stars at approximately equal temperatures and thus 
we expect the energy deposition of eqn. (3), which scales like oc T 9 , to be roughly constant 
for stars differing in mass by 5 %, which is a reasonable range for observed neutron star 
binaries. 

In conclusion, we find that inter-star neutrino annihilation is a significant component 
of the total neutrino annihilation in a neutron star binary system and will thus be a major 
contributor to the available energy of coallescing binary systems. We provide formulas 
for calculating this energy deposition. Future work will use 3-D relativistic hydrodynamic 
calculations to take this effect into account and thus improve the calculations of Salmonson 
et al. (2001). In particular the inter-star annihilation will lead to non-isotropic outflow. 

This work was performed under the auspices of the U.S. Department of Energy by 
University of California Lawrence Livermore National Laboratory under contract W-7405- 
ENG-48. 
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Fig. 1. — Diagram representing the geometry and defining quantities used in Eqns (9, 11, 
12). For clarity, this illustration shows only the special case where the neutrino trajectories 
are coplanar (i.e. 4>i — 02 = 0). 
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Fig. 2. — A contour plot of a slice through the stars of the angular dependence of the 
inter-star neutrino energy deposition density, Ointer, given by eqn. (9) with Eqns (11, 12), 
with stellar radii 10 km and separation 30 km. There are ten 1.1 increment contours starting 
from black, (0,1.1), and increasing to white. The neutron stars are represented by grey disks. 
This figure is meant to illustrate the function O^ter and thus does not include eclipsing of 
neutrino paths by the neutron stars which would attenuate this function around the far sides 
of the neutron stars. Since the effect is already small in these regions, this modification is 
negligible. 
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Fig. 3. — A plot of the region between the two neutron stars along the line connecting the 
two star centers. The left (right) edge of the plot is the surface of the left (right) star. 
The solid line is the inter-star deposition Oi nte r given by Eqn (9, 11, 12) and shown also in 
Figure 2. The dashed lines are the single star angular deposition functions S mgic (Eqn. 13, 
14, Newtonian; M — 0) for each star of radius 10 km and separation 30 km. Note that 
©single decays rapidly with distance from the surface of each neutron star while O^ter is more 
gradually varying between the stars. 
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Fig. 4. — Ratio of integrals over space of inter-star energy angle factor (Eqns 9, 11, 12) to 
single star angle factor (Eqns 13, 14, Newtonian; M — 0). This is evaluated for star radii 
10 km over a range of half-distances, d 2 , between the star centers. The turn-down for small 
separations (ok < 12 km) is due to the paucity of volume between the stars. 
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Fig. 5. — A ratio of the entire energy deposition due to inter-star neutrino annihilations, 
Qinter, to that from single star annihilations, Qsmgie, is shown for isotropic stellar radii R = 7.5 
km for a range of stellar half-separations, d,2, and for several stellar masses M. M = 
corresponds to the Newtonian, zero-gravity case. Note that for M = 2 km 1.4M and 
d 2 = 10 km inter-star energy deposition, Qi nter , equals single star energy deposition Qsingic- 
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Fig. 6. — A comparison of the ratio QcTi/QNewt of general relativistic (eqn. 26) to Newtonian 
(eqn. 15) energy deposition. This ratio is shown for inter-star Qmter and single star Qsmgie 
energy depostions for a range of stellar masses and with isotropic stellar radius R = 7.5 
km and half-separation di = 15 km. As mass increases, both components are augmented 
by relativistic effects. For reference, the energy deposition is shown for an isolated neutron 
star, Q one star , i.e. the same as Q S ingio, but with the companion star mass set to zero. The 
depositions in a binary system (Qinter, Qsingic) are augmented over Q onc star because of the 
deeper, two-star potential well. 



